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Abstract 



Two types of pseudoscalar quarkonium electromagnetic decay processes, i.e. decay to a lepton 
pair, and to a lepton pair plus a photon (Dalitz decay), are analyzed at the leading order in 
NRQCD expansion. The former type of processes, highly suppressed in the Standard Model, have 



> 
(N 
oo 

■ been hoped to act as the sensitive probes of the possible new physics. The latter type of processes 



generally possess much greater decay rates than the former, owing to several conspiring factors. 
The recently launched BES-III program, with 10^ ??c samples to be anticipated in the coming years, 
may be able to observe the Dalitz decays r]c — > e^e~7 and r/c — > fj,'^ which have branching 
ratios of order 10~^. When the radiated photon becomes very soft, the Dalitz decay events will 
be experimentally tagged as the exclusive lepton pair events. It is found that, those quasi-two- 
body events that arise from rjc — > e~^e~'y with photon energy less than the minimum sensitivity 
of the electromagnetic calorimeter, can vastly outnumber the literal r]c e^e~ events, however 
this amplification is still not dramatic enough for the BES-III experiment to establish these events. 
Consequently, the expectation of looking for new physics signature in the rjc l^l^ channel is 
obscured, unless the contamination from rfc I'^l^j has been taken into account carefully. 

PACS numbers: 12.38. -t, 12.38. Bx, 13.20.-v, 13.20. Gd, 13.40. Hq 
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I. INTRODUCTION 



Charge-neutral pseudoscalar meson decay to a lepton pair, P ^ l^l , has long been 
an interesting topic. In particular, for P to be a light pseudoscalar meson (vr^, rj, ■)]'■■■), 



extensive theoretical and experimenta 



study of 77° — ^ e^e in 1959 



efforts have been conducted since Drell initiated the 



ii,i2i,y,y,B,y,0,y,3. 



Studies of these decays can offer insights 



into the nonperturbative structure of the pseudoscalar mesons, in particular, help one to 
glean more knowledge about the P7*7* transition form factor. This type of electromagnetic 
decay processes are suppressed by two additional powers of a with respect to P ^ 77, 
and also penalized by helicity conservation. As a consequence, the decay probabilities are 
generally very tiny within the Standard Model (SM), rendering experimental detection rather 
challenging. On the other hand, the rareness of these decay processes may turn into a virtue, 
that is, they might be utilized as the sensitive probes of possible new interactions beyond 
SM. 

Another intimately related type of electromagnetic decays are P —>■ These kinds 

of pseudoscalar meson decay processes are of use to extract the information about the P7*7 
form factor. Despite radiating off an extra photon, these bremsstrahlung leptonic decays 
in general occur much more copiously than P l^l^ . For example, the so-called Dalitz 



decay process 



10|, TT*^ — e^e 7, which has been observed decades ago, has a branching 



ratio of (1.198 ± 0.032)% ll|]. This is more than five orders of magnitude greater than the 



branching fraction of vr e^e , (7.48 ± 0.38) x 10 , which was recently measured in the 



KTeV E799-II experiment at Fermilab [12]. This striking disparity can be attributed to 
a number of facts, that such bremsstrahlung leptonic decay processes are suppressed with 
respect to P ^ 77 by only one additional power of a, suffer no helicity suppression, and 
also enjoy coUinear enhancement brought in by photon fragmentation to a lepton pair. 

In this work, we aim to investigate both types of leptonic decay processes for P to be 
a pseudoscalar heavy quarkonium state, i.e., rjQ —>■ l^l" and rjQ —>■ l'^l~'-f {Q = c,b) ^. 
One strong incentive stems from the experimental side. For instance, the recently launched 
BES-III experiment, plans to accumulate an unprecedentedly large data set of charmonia, 
e.g., 10^^ J/i/j and 3 x 10^ ^'(25') in the coming years |l^. An enormous number of 7]c and 



^ In literature, the term Dalitz decay has been specifically reserved for tt'^ e+e 7 [10]. In this work we 
generalize the use of this term, i.e. we also use it to refer to any of the rjQ l^l^j process. 
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ric{2S) are expected to be produced via the radiative transitions from these J/ip or 4'(2S) 
samples. Analogously, the scheduled Super Belle experiment, will also be capable of col- 
lecting a tremendous number of rjh samples. Further, an even larger data set of pseudoscalar 
quarkonia are expected to be produced at the CERN Large Hadron Collider (LHC), though 
the copious backgrounds in hadron machine renders the detection of such rare decays rather 
challenging. In any event, it seems, at the current time, not of only academic interests 
to assess the observation potentials of these rare electromagnetic decay processes in the 
forthcoming experiments. 

From a theoretical perspective, it is also worthwhile to study these rare electromagnetic 
decays of pseudoscalar quarkonia, to enrich our knowledge about heavy quark physics. A 
heavy quarkonium state, being a heavy-quark heavy-antiquark pair tightly bound via the 
strong interaction, is the best understood among all types of hadrons. In sharp contrast to 
light meson decay, which must be analyzed by some nonperturbative tools, quarkonium decay 
can be accommodated in the perturbative QCD framework, owing to the condition itlq 3> 
^QCD- Indeed, annihilation decays of heavy quarkonium, especially the electromagnetic 
ones that we plan to investigate, can be systematically tackled by the modern effective- 
field-theory formalism, the nonrelativistic QCD (NRQCD) factorization approach [14? ]. 
In a quarkonium electromagnetic decay process, this factorization approach allows one to 
systematically separate the hard quantum fluctuation of order heavy quark mass m from 
the low-energy contributions of order mv or smaller, where v signifies the typical velocity of 
Q or Q in a quarkonium. Empirically, ~ 0.3 for charmonium, and 0.1 for bottomonium. 

We note that, both types of leptonic decay processes of r/g have already been p artly 
investigated by different authors. For instance, rjq has been studied in Ref. jl, 15, IGj . 



and the Dalitz decay rjq — > was considered in 17|, [iSj. Nevertheless, a comprehensive 

analysis based on the NRQCD factorization approach is still lacking ^. For this reason, we 
feel that it might be rewarding to revisit these two processes from this angle. We will work 
at the leading order in NRQCD expansion only. However, the systematics of this approach 
renders future implementation of higher-order corrections possible. 



^ Ref. 4] studies the rjQ l^l^ process within a bound state quark model, which nevertheless may be 
viewed as a primitive version of the NRQCD approach. Ref. [l^ employs the NRQCD factorization 
explicitly, but only studies a few processes using numerical recipe. Both of works completely neglect the 
weak interaction contribution, which turns out to be inadequate for ry^ decay. 
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We summarize the main outcome of this work. We confirm the analytic expression for the 
electromagnetic contribution to the decay rjQ l^l^ , which was first reported in Ref. 4|. 
We also include the often-omitted weak interaction contribution arising from Z° exchange. 
It is found that, although the weak interaction plays a negligible role in the leptonic decay 
of ?7c, its effect can become important in rjt, decay, especially for r^t t^t^ ^. Despite this, 
the net SM predictions to the decay rates are still too suppressed for these processes to be 
observed experimentally in the foreseeable future. 

We also perform a comprehensive study of numerous pseudoscalar quarkonium Dalitz 
decay processes, and verify that such decays generally possess much more enhanced decay 
probability than t]q — * l^l^. In particular, it is found that rjc e~^e~'y, fi'^fi~'~f, with 
branching fractions of order 10~^, might have bright prospect to be observed at the BES- 
III experiment. We also analyze the energy distributions of lepton and photon in these 
Dalitz decay processes. We hope future measurements of these energy spectra can test our 
predictions critically. 

We emphasize that any realistic electromagnetic calorimeter is limited by the finite sen- 
sitivity to detect soft photons. This indicates that, the Dalitz decay event t]q l^l^'-f will 
fake the literal tjq l^l^ event, if the emitted photon is too soft to be registered by the elec- 
tromagnetic calorimeter. Therefore, the hope of seeking new interaction beyond SM in the 
rjQ l^l^ channel becomes obscured, unless the contamination from the respective Dalitz 
decay is thoroughly understood and incorporated in the analysis. Our study shows that, the 
number of quasi-two-body events that arise from rjc e"'"e~7 with photon energy restricted 
to be less than 20 MeV, can easily surpass that from rjc — > e"'"e^ by two orders of magnitude, 
but such enhancement is still not significant enough to warrant the establishment of such 
events at BESIII experiment. 

The remainder of the paper is distributed as follows. In Section [Tll we calculate the 
pseudoscalar quarkonia decays to a lepton pair at the leading order in NRQCD expansion, 
including both contributions from QED and weak interaction. We also compare our results 
with the previous ones. In Section IIIIl we present a detailed calculation of the pseudoscalar 
quarkonium Dalitz decay processes at the leading order in NRQCD expansion. The inclusive 
energy spectra of the photon and the lepton are presented, and the QED fragmentation 



The weak interaction effect for rjc I I has also been considered in [16|. 
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function for a photon to split into a lepton can be extracted. We also obtain a succinct 
expression for the decay branching fraction integrated over the full three-body phase space. 
In Section \IV\ a comprehensive numerical analysis for numerous processes of pseudoscalar 
quarkonia decays to a lepton pair and the Dalitz decays are made. Particular attention 
is paid to the interference pattern between QED and the Z°-exchange contributions to 
the former type of processes. We assess the observation prospects of both types of decay 
processes. We also investigate to which extent the Dalitz decay events will fake the exclusive 
lepton pair events. In Section |Vl we summarize and present an outlook. 

II. PSEUDOSCALAR QUARKONIUM DECAY TO A LEPTON PAIR 




al) a2) b) 



FIG. 1: Lowest-order diagrams for the processes rjq l^l^ [a] and rjQ \b]. In the former 

process, we represent the QED contribution by al), the weak interaction contribution by a2). The 
crossed diagrams for al) and h) have been suppressed. 

The purpose of this section is to derive the analytic expressions for the amplitude of 
rjQ in the lowest-order NRQCD expansion, by taking only the SM interactions into 

consideration. That is, we will consider the contributions from both electromagnetic and 
weak interactions, and investigate their interference pattern. 

A. The electromagnetic contribution to r]Q —>■ 

In this subsection we focus on the QED contribution to this process. Unlike the leptonic 
decay of vector quarkonium such as J/V', a pseudoscalar quarkonium cannot directly decay 
to a lepton pair through annihilation into a virtual photon at tree level. At the lowest order 
in electromagnetic and strong couplings, this process proceeds through the one-loop QED 
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box diagram, as depicted in Fig. [T^) 

The validity of the NRQCD approach rests upon one of the key characteristics of quarko- 
nium, that both of its constitutes move non-relativistically in the quarkonium rest frame, 
so the quark relative velocity can serve as a small expansion parameter of the theory. For 
this reason, NRQCD would be a very poor framework to describe light mesons such as vr, 77, 
77', etc. Nowadays the NRQCD approach has been accepted as the standard tool to analyze 
heavy quarkonium decay and production processes. Since hard reactions involving quarko- 
nium necessarily probe the scale of order heavy quark mass, by appealing to the asymptotic 
freedom of QCD, NRQCD approach allows one to put the amplitude in a factorized form, 
i.e., the sum of the products of perturbatively calculable short- distance coefficients and 
nonperturbative but universal NRQCD matrix elements. 

The calculation of quarkonium decay at the LO in NRQCD expansion is standard. We 
assume the quarkonium state composed of a heavy quark Q and its antiquark Q, and abbre- 
viate the QQi^S^^) state by t]q. We assign the momenta carried by r^g, /~, as P, pi, p2- 
At the LO in velocity expansion, one can routinely obtain the amplitude for r/g by 
first computing the amplitude for Q{^)Q{^) — > l~ {pi)l^ {P2) ■, enforcing Q and Q to carry 
equal momentum, then projecting it onto the intended ^S'^'* state. 

There are totally two lowest-order QED diagrams for this process, one of which is shown 
in Fig. [T]-al). The other undrawn diagram can be obtained from it by reversing the fermionic 
arrow in either the quark or the lepton line. By C-invariance, both of diagrams yield the 
identical results. After some straightforward calculation, we can express the electromagnetic 
decay amplitude at the LO in NRQCD expansion as 



where Nc = ?> i?> the number of colors, e; and eg signify the electric charges of charged 
lepton and heavy quark in units of |e| (e; = —1, Cc = 2/3 and Ch = —1/3), and a is the fine 
structure constant. The nonperturbative factor ?/'^^(0), is the wave function at the origin 
for the riQ state, which can be identified with the LO NRQCD matrix element, m; and mg 
denote the masses of lepton and quark, respectively. At the current level of accuracy, it is 

^ The same annihilation type of diagram, with the lepton replaced by the light quark and the photon 
replaced by gluon, has been hypothesized to account for the rjc — ?7(?7') mixing [191] . 




(1) 
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legitimate to treat mg and Mrf^jl interchangeably. The dimensionless function / encodes 
the effect of loop contribution, and is normalized in such a way that it depends on the ratio 
of lepton mass to quark mass at most logarithmically. 

Some remarks on the traits of this process are in order. Because the lepton pair must 
form a ^5*0 state to conserve angular momentum, a 75 is expected to be sandwiched between 
the leptonic spinors in the decay amplitude. Indeed, with the aid of Dirac equation, one 
easily verifies that the bispinor ^(^1)75^(^2) exhausts all the possible Lorentz structures. It 
is well known that leptonic decays of pseudoscalar meson suffer from the so-called helicity 
suppression, and would be strictly forbidden when lepton mass set to zero, hence there 
should be an explicit factor of lepton mass appearing in the amplitude ^. Note that the 
/ function diverges with mi only logarithmically, so equation ([T]) is compatible with the 
helicity suppression mechanism. 

It turns out that the evaluation of the box diagram can be reduced to evaluating a three- 
point one-loop integral. The encountered loop integrals are both ultraviolet and infrared 
finite, hence may be directly computed at four spacetime dimension. After some efforts, we 
can obtain the closed form for the / function: 



4 \1-I3 \1-I3) 12 n l+lV 2 \l-l3 



(2) 



where Li2 is the dilogarithm. We have introduced r = -^tt- = —3-, and /5 = v 1 — r is the 

M^^ nig 

velocity of the outgoing lepton in the rig rest frame. The shape of this function is shown 



once the 



in Fig. [2j Note Eq. ([2]) is identical to the R function given in equation (12) of 
Spence function $(x) there is identified with — Li2(— x). It may be worth pointing out that, 
in order to reach the compact expression given in ([2]), we have made use of the following 
somewhat inapparent relation ^: 



^ We stress this phenomenon not only pertains to pseudoscalar meson leptonic decay. The decay Xco ~^ e~^e~ 

should also be forbidden in the zero-We limit, again due to the conflict between angular momentum 

conservation and helicity conservation in massless QED. 
^ This relation can be proven with the help of the identity Li2 ( f ) = Li2 (x) + Li2 ( ^ ) + Li2 ( 2{x-i) ) ~^ ^^"^ ^•'^ ~ 

1) + ^ ln^(2— 2a;). Substituting x — I — (3 and x — —jj^ into this identity separately, taking the respective 

difference, one then obtains the desired answer. 
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FIG. 2: The profile of the function f{r). Sohd lines represent the exact results given in (l2|), and 
the dashed ones represent the asymptotic expressions as given in (f3a|) and (i3bl) . 

The occurrence of the imaginary part in ([2]) is linked with the "unitarity bound" of the 
branching ratio, which is obtaining from tieing the amplitudes of rjq — > 77 and 77 
together according to the cutting rule. 

In Nature heavy quarks are generally much heavier than leptons, especially for the first 
two generations. It is then useful to know the asymptotic behavior of the / function in the 
r — i> limit: 



As anticipated, /(r) depends on r only logarithmically. From Fig. [21 one can see that the 
function fasym{r) already constitutes a rather good approximation as r < 0.1. To reproduce 
these asymptotic behaviors more efficiently, one may appeal to the method of region {2^ 
by dismembering the original loop integration into the sum of integrations from different 
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regions, e.g., hard, soft, collinear-to-/~ and collinear-to-/"'" in our case j2l|]. One may identify 
the double logarithm in fl3a|) that originates from the overlap between coUinear and soft 
singularities. 

In passing we remark on one peculiarity reported in a recent calculation of the same 

I 

process that employing the light-cone approach ^J!]. In that approach, the amplitude is 
expressed as the convolution of a hard-scattering part with the light-cone distribution ampli- 
tude of rjQ. The hard part there is found to scale as lnr/^/r in the limit r ^ 0. This infrared 
sensitivity is theoretically disastrous because it diametrically conflicts with the requirement 
of helicity suppression- that this process should be strictly forbidden for a massless lepton. 
This nuisance, if confirmed to persist, may indicate that, the light-cone operator product 
expansion (OPE), which underlies the calculational framework in [iqI, may no longer be 
suited to describe heavy quarkonium decay. The light-cone expansion is usually formulated 
as an expansion in terms of a small energy scale like quark mass over a large momentum 
transfer. It should be appropriate for a high-energy exclusive heavy quarkonium production 
processes, as the large momentum transfer scale can be identified with the center-of-mass 
energy of the reaction, which may indeed be much greater than the heavy quark mass. 
However in the process at hand, mg itself already acts as the highest energy scale, so it is 
difficult to imagine the actual meaning of the light-cone expansion here. In a space-time 
picture, the slowly-moving heavy quark and antiquark typically experience an instantaneous 
strong force, so they are typically separated by a distance of order l/rnqv but local in time, 
which is quite far from a light-like separation. In our opinion, NRQCD factorization ap- 
proach, which is closely related to a local OPE by treating 1 / mq as an expansion parameter, 
provides the most natural and economic framework to account for the heavy quarkonium 
decay, in particular for the processes considered in this work. 

It is also interesting to look at the alternative limit r ^ 1 (/3 — 0), where the mass of 
rjQ is just sitting at the threshold of twice lepton mass. From ([2]), one finds the following 
limiting value for /: 

/(I) = -2 + 21n2-i7r, (4) 



which can also be seen in Fig. O As noted in Ref. this is a well-known result, which is 
responsible to the two-photon annihilation contribution to the hyperfine splitting between 



the orthopositronium and parapositronium 



221 ] (see also 



23 
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B. The weak- interact ion contribution to rjQ ^ I 

Neutral quarkonium decay is normally not a good place to look for the trace of weak 
interaction, which is generally overshadowed by the strong and electromagnetic interactions. 
However, for the processes at hand, the LO QED amplitude has to proceed at one loop order, 
but the weak interaction contribution instead can start at tree level, so it is not inconceivable 
that the weak interaction may play some role for some of these processes. 

There is only one s-channel diagram that contributes to ?7q from Z° exchange, as 

depicted in Fig. [T]a2). Only the axial vector coupling of Z^ff contributes to this process. 
At the lowest order in the velocity expansion, the weak amplitude can be expressed as 

^WEAK = 2v/2iV^^ MPlhMP2)], (5) 

y/fn^ sm Ow COS'' Ow 

where 9w is the Weinberg angle, and 5^^, g'^ denote the weak axial charge of charged lepton 
and quark, respectively. The weak axial charge of a fermion is equal to its 3rd component 
of weak isospin ((7^ = — ^ for I = e,/i,r; g'^ = |, fl'^ = ~|)- We have neglected the rjQ 
mass as well as the width of in the Z^ propagator, since they are much smaller than Z° 
mass for Q = c, b. Therefore the Z° exchange can be effectively mimicked by a four-fermion 
contact interaction. Note there is an explicit factor of mi in Eq. ([5]), again due to helicity 
suppression, very similar to what occurs to 71^ l^vi- 

Comparing ^ with ([T]), one clearly sees that, though the weak interaction amplitude is 
suppressed by a factor of mg/Mf relative to the electromagnetic one, it suffers less sup- 
pression by one power of a. For r^^ decay to a lepton pair, these two competing effects may 
become comparable in magnitude. As a consequence, in order to make reliable predictions, 
the weak interaction effect becomes indispensable and must be included. 

C. The Standard Model prediction to ryg — > / Z+ 

Substituting Eq. ([T]) and Eq. ^ into 

^^^[riQ ^ = ^ I^EM + ^^weakI' , (6) 

we then get the desired partial width expected in the frame of SM. 
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It is convenient to introduce the normalized decay rate of r^g — / defined as Eq. 
normalized to the partial width of rjq — * 77: 

T^M^rjQ ^ 1+1-] 



To 



a 

2^ 



pr 



fir)-g'i 



le2 a2 



(7) 



to 



where we have substituted ei = —1, and used the relation Mz = ( ^'^ V^"^ . „ ^ — 5—, 
condense the expression for the scaled weak amplitude {Gp is the Fermi coupling constant). 
An important feature is that the relative importance of weak interaction contribution in- 
creases with Mq. The partial width of rjq 77 is given at the lowest order in v and 
as- 



To = ^[vq^ 77] 



m 



Q 



(8) 



The advantage of introducing the R ratio in ([7]) is that the nonperturbative factor ipr^d^) 
cancels out in the ratio, and some portions of QCD radiative and relativistic corrections 
to both rjQ l^l^ and rjQ 77 may largely cancel. Therefore, for a LO calculation like 
this work, using the R ratio instead of the branching ratio is presumably more appropriate. 
In the situation where the r]Q di-photon decay has been experimentally measured, one can 
directly obtain the B[rjQ /+/^] by multiplying BcxplVQ ~^ 77] with the predicted R ratio. 



III. DALITZ DECAYS r/g ^ l+l'-f 

It is an experimental fact that, for light pseudoscalar mesons, the branching fractions 
of the Dalitz decay processes are several orders of magnitude greater than those of the 
respective leptonic decays. Certainly, it is worthwhile to examine whether the same pattern 
also holds for pseudoscalar quarkonium decays or not. 

We note that, an analogous pseudoscalar quarkonium strong decay process, i.e. rjQ qqg, 
has been previously studied by several authors by retaining a nonzero mass for q. The 



analytic expressions for the energy distri 
decay rate, has been presented in Ref. ] 



3utions of the gluon, quark and the integrated 
251 ]. In the following, we will independently derive 
the corresponding energy spectra of the photon, lepton and the integrated decay rate for the 
rjQ Dalitz decay process. When inserting a proper color factor, the exact agreement is found 

11 



15 



between our Eqs. (fT3l) . (fT^ . (fT6|) and Eqs. (4),_ (5), (6) in Ref. [25l]. The process rji, ccg 
has also been studied numerically in 



26 



27| 



A. Squared amplitude of quarkonium Dalitz decay 

In contrast to the rare decay r]Q l^l^, the Dalitz decay tjq — > can start at 

tree level in QED ^. There are in total two diagrams, one of which is depicted in Fig. [Tf)). 
The momenta of r/g, l" , 7 are assigned as P, pi, p2, ps, respectively. This process first 
proceeds through 7]q —>■ 77*, and the virtual photon then fragments into a lepton pair. If the 
lepton is much lighter than the quark, one expects that the decay rate is dominated by the 
kinematic configuration where the invariant mass of lepton pair is close to its minimum, 2m/. 
In other word, it is the nonzero lepton mass that cuts off the potential coUinear singularity. 

It might be worrisome that the amplitude may diverge in another kinematic configuration, 
i.e., where the photon becomes very soft, and the lepton and anti-lepton move nearly back- 
to-back with equal momentum. It is well known that in J/ip decay to l^l'j, infrared 
divergence does arise in the long wavelength limit of photon. This divergence is in turn 
canceled by including the virtual correction to J ftp —>■ guaranteed by the Bloch- 

Nordsieck theorem ^. However, such mechanism obviously does not apply to our case to 
sweep the potential infrared divergence. 

A little thought reveals that Fig. ^) must be regular in the — >■ limit. It is most 
transparent to see this in the context of nonrelativistic effective theory. To describe an 
almost on-shell quark interacting with a soft photon, one is justified to use nonrelativistic 
QED (NRQED). After the incoming Q emits a soft photon, it has to annihilate with Q into 
a virtual photon that splits into a lepton pair. Therefore, the emission of this soft photon 
has to flip the spin of Q, i.e., effectively induces a magnetic dipole transition, to convert the 
QQ pair from the initial ^5*0 state to a ^Si state. The effect of this soft-photon emission is 
accounted for by the operator ^^^^cr ■ B'^^tp, where ip denotes the Pauli spinor field for Q 
in NRQED. The occurrence of the magnetic field strength, S*^™, which brings forth a factor 



^ Weak interaction contribution can also start at tree level, nonetheless is completely negligible for both 
rjc and r/b Dalitz decays. 

® Experimentally, J/ip — s- events are selected by requiring the photon energy is greater than the 



characteristic detector sensitivity, say, 100 MeV in the Fermilab E760 experiment [28|. Those three-body 



decays with photon energy less than 100 MeV are tagged as the J/ip l^l events. 
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of ps, will protect against the infrared singularity arising in the quark propagator, therefore 
the amplitude is infrared finite ^. 

At the lowest order in strong coupling and in t>, the decay amplitude can be routinely 
obtained: 

where e-^ signifies the photon polarization vector. We have taken the undrawn diagram into 
account, which yields identical contribution as Fig. ^) owing to the C-invariance. 

From ^ we are reassured that, the numerator of the amplitude does contain a factor of 
photon momentum, which is crucial to tame the infrared singularity arising from the quark 
propagator. It is also worth noting that, since the lepton and anti-lepton directly come from 
the photon fragmentation, they necessarily form a ^Si state. There is no helicity suppression 
mechanism affiliated with this process, consequently no factor of mi manifests in ([9]). 

For notational abbreviation, we introduce three dimensionless energy variables as Xi = 
2Ei/M^Q = 2P ■ p,/P^ {i = 1, 2, 3), where = M^^, = p° is the energy of each final 
state particle in the rjQ rest frame. We also adopt the same definition as in the previous 
section, r = imf/M^^. 

It is straightforward to square the amplitude and sum over the polarizations of final-state 
particles. We then obtain 



- xi) + 0:2(1 - X2)] xs - 2(1 - a;i)(l - X2) + \x\ 



xl{\-x,Y ' 

which is symmetric under the interchange between x\ and xi- As expected, each of the three 
terms in the numerator scales as x\ in the infrared limit X3 — 0. 

Energy conservation demands that xi + X2 + 0:3 = 2. Any kinematic invariants can be 
expressed in terms of any pair of these three scaled energy variables. In the following we 
shall derive the inclusive energy distributions of /~ and 7 separately, hence it is convenient 



^ By the similar reasoning, it is easy to see that the decay xqj ^ l^l^j does develop an infrared singularity 
in the soft-photon limit, since the corresponding NRQED operator governing the soft-photon emission is 
of the electric dipole type, i.e. ~2mQ V'^ [-P ' ^'^"^ + ' -DjV': and there appears no factor of photon 
momentum to kill the corresponding one in the denominator. 
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to choose Xi and xs as the two independent variables in expressing the differential three- 
body phase space. The dimensionless energy variable of X2, has been eliminated in 
the amplitude squared. The energy spectra of 7 and /~ can be obtained by choosing the 
different order of two-fold integrations over Xi and X3. As a consequence, the decay rate can 
be expressed in the following two different ways: 

/.l-r I-X+ 

r[r^Q^r/+7] = dxs J dx,J2\^\' (11a) 

/•I rx^ 

/ dxi / dx3j2\-^\'^- (11^) 



32(27r)^ 



The integration boundaries for the outer-layer integrals have been labeled explicitly; for the 
inner-layer integral, the upper and lower boundaries can also be readily inferred: 



^ 2-3:3 , jl-r-xs 



2 — xi =F V ^1 ~ ^ 
B. Inclusive photon energy spectrum 

It is straightforward to deduce the energy distribution of the photon by integrating over 
the variable Xi in (11 lap : 

dRjrjQ ^ 7(0:3) +X] _ 1 dT[r]Q ^ 7(3:3) + X] 
dx'i To dx-i 



_ a 3:3(2 + r - 2a:3)v/l - r - X3 , , 

- (1-3:3)5/2 ■ ^ 

As in dealing with rjq decay to a lepton pair, we introduce the scaled energy distribution of 
7, dR{x3)/dx3, the differential decay rate of rjq ^ •y + X normalized with respect to Fq, the 
partial width of riQ 77 given in ([8]). This distribution clearly vanishes at both the lower 
and upper ends of the scaled photon energy, i.e. X3 = and X3 = 1 — r. For r <^ 1, the 
spectrum is generally featureless and negligible in most of the region, except a sharp peak 



rises near the upper end of X3, which centers at 0:3 ~ 1 — ^121+1^ ^ 1 — 1.40r, and the peak 
height ~ O.lla/r. This clearly indicates that, the decay rate is dominated by the kinematic 
configuration where the outgoing lepton pair carries a small invariant mass, owing to the 
collinear enhancement. 
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There is also motivation to inspect the lower end of the photon spectrum more closely. 
It is an experimental fact that a realistic electromagnetic calorimeter can detect photons 
only down to some minimum limiting energy E^^ut- If the photon becomes very soft, it 
will not be properly registered by the electromagnetic calorimeter. In this respect, those 
t]q l^l^'j events with a very soft photon will, from the experimental perspective, mimic 
the respective tjq event. We learn from (fT3|l that, at small X3, this differential decay 

rate becomes enormously suppressed relative to that in the 0:3 —> 1 limit, and scales linearly 
with X3, ~ — r(2 + r)x3. Imagine we impose a realistic cutoff x^cut on the photon 

energy. Those three-body decay events will be correctly recorded as Dalitz decay events 
when the fractional photon energy is greater than x^cuu in comparison, those three-body 
events will be tagged as the lepton pair events when X3 < x^cut- The R value integrating over 
x^ from to this cutoff will be approximately for small r. Although this is a tiny 

fraction, the chance exists that for some Dalitz decays, it might still be much greater than 
the extremely small decay ratio of rjQ l^l^ as given in (I7j). If this is the case, what are 
experimentally recorded as the rjQ l^l^ events in fact receive the bulk of contributions 
from the three-body Dalitz decays. As we shall see, r]c — >■ e~e"'" constitutes such a very 
example. 



C. Inclusive lepton energy spectrum and photon-to-lepton fragmentation function 

We are also interested in the inclusive energy distribution of the lepton. This can be 
obtained by integrating over the variable 0:3 in (lllbl) : 

dR[r]Q ^ l-{xi)+X] _ 1 dT[r]Q ^ /"(xi) +X] 




(14) 

As in Eq. f fT3|) . we also define a dimensionless distribution, dR{xi) / dxi, the differential decay 
rate of t]q l~{xi) + X normalized with respect to Fq. It is easy to see that this energy 
spectrum vanishes at both end points: Xi = ^/r and 1. In sharp contrast to the energy 
spectrum of 7, the spectrum of /~ is more evenly populated in the whole region of Xi. In 
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FIG. 3: Normalized energy distributions of leptons in numerous i]c{r]'^) and r]h Dalitz decays. The 
energy spectrum of r in t]'^ Dalitz decay, populating a rather narrow region near xi = 1, is very 
much suppressed with respect to the other ones, so is not displayed in this figure. 

Fig. [3l we show the energy distributions of different species of /~ in various Dalitz decays of 
Tie, ?7c and r]f, 

As r <^ 1, the formalism of fragmentation function can adequately account for the lep- 
tonic energy distribution. One can readily identify the fragmentation function D^f^i~{z), by 
expanding (HM in powers of r, retaining only those nonvanishing terms in the limit r ^ 0, 
and dividing them by 2 to compensate the fact that each photon in rjQ 77 can fragment: 



(15) 



As expected, the coefficient of the collinear logarithm is nothing but the Altarelli-Parisi 
splitting kernel P^^i-. We have checked that, for each Dalitz decay channel tjq l^l^'j 
(except I = t), this fragmentation function approximates the exact spectra f|T4l) quite well, 
and also gives a satisfactory account of the corresponding integrated decay rates. 



^'^ Note that the quark energy distributions in various r/g qqg processes, have not been correctly displayed 
in Fig. 3 of due to some input mistakes. 
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D. The integrated decay rate of rjQ — > l^l 7 

It is also desirable to know the total decay rate of t]q by integrating the 

amplitude squared in flTUl) over the full three-body phase space. The integrated decay ratio 
can be readily deduced by starting from either the photon spectrum in (1131) or the lepton 
distribution in (IT^ . The final answer admits a particularly succinct form: 



(16) 



Its limiting behavior in small r reads 

[integrated] = ^ jln + 2 In 2 - ^ + C(r) | . (17) 

As expected, the greater the disparity between quark mass and lepton mass is, the larger the 
integrated R value becomes due to the increasing logarithmic enhancement. This asymptotic 
expression can also be easily reproduced by integrating twice of the fragmentation function 
D^f^i~{xi) in ( |T5|) over the entire range of the energy fraction Xi of /~. 

It is also instructive to look at the alternative limit r ^ 1, where the r/g is barely heavy 
enough to disintegrate into a leptons pair at rest plus a zero-energy photon: 

i?[integrated] = ^ + 0(/3^)} , (18) 



where (3 = y/1 — r. This somewhat unexpectedly severe suppression, will be useful for us to 
understand the very small decay rate of rjc{2S) r^r^-y. 

IV. PHENOMENOLOGY 

In this section we will explore the consequences of the formulas presented in previous 
sections. In particular, we will analyze numerous leptonic and Dalitz decay channels of 
pseudoscalar charmonia and bottomonia, to assess their observation potentials in the current 
and forthcoming high-energy collision facilities. 

In Table [H we tabulate numerous predictions for rjc, ric{2S) and rji, decays to all possible 
species of lepton pairs. When evaluating the mass ratio r, we take the precisely known lepton 



masses 



ll|: TUe = 0.511 MeV, = 105.66 MeV, and = 1776.84 MeV, and use the 



central values of the less precisely measured pseudoscalar quarkonium masses [ll|: M, 
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TABLE I: The values of r, the scaled amplitudes for electromagnetic and weak interactions, and 
the respective R ratios for various pseudoscalar quarkonium decays to a lepton pair. is the 

R ratio that includes both QED and weak contribution, as indicated in d?]); R^^ denotes the 
corresponding R ratio by retaining only the QED contribution. For some of the decay channels. 



NRQCD predictions to R^^ have also been given m 



Decay modes 


r 




/( 


r) 


Q V2GpM^^ 
w 










77c e+e" 


1.18 X 10" 


-7 


58.67 - 


27.24i 


-0.39 


1.33 X 10" 


-9 


1.31 X 10" 


-9 


ric /^^^~ 


5.03 X 10- 


-3 


5.30 - 


10.51i 


-0.39 


1.88 X 10" 


-6 


1.82 X 10" 


-6 


7],{2S) ^ e+e- 


7.90 X 10" 


-8 


61.76 - 


27.87i 


-0.58 


9.79 X 10" 


10 


9.63 X 10" 


10 


7]c{2S) n+fi' 


3.38 X 10" 


-3 


6.27- 


11.13i 


-0.58 


1.49 X 10" 


-6 


1.42 X 10" 


-6 


ri,{2S) ^ T+r- 


0.955 




-0.61 - 


- 3.19i 


-0.58 


5.78 X 10" 


-6 


6.35 X 10" 


-6 


rib e+e" 


1.19 X 10" 


-^8 


77.59 - 


30.85i 


15.35 


2.23 X 10" 


10 


3.07 X 10" 


10 


% 


5.07 X 10" 


-4 


11.98 - 


ll.lOi 


15.35 


4.68 X 10" 


-7 


1.29 X 10" 


-6 




0.143 




0.19 - 


5.52« 


15.35 


1.09 X 10" 


-5 


9.74 X 10" 


-5 



2980.3 MeV, M^^(25) = 3637.0 MeV, and M^, = 9388.9 MeV |29|. For the electromagnetic 
and weak couplings, we take a = 1/137 (for simplicity, we have neglected the running effect 
of fine structure constant), and Fermi coupling constant Gp = 1.166 x 10"^ GeV"^. We 
have taken the electric charges of quarks to be Cc = 2/3 and Cb = —1/3, the weak axial 
charges of quarks to be g% = 1/2, g\ = —1/2, respectively. 

Listed in Table [T] are various normalized R values. Since each leptonic decay of t]q is 
severely suppressed with respect to rjQ 77, i? -C 1 certainly is expected. If one wishes 
to convert these R ratios to the corresponding branching fractions, one should multiply 
them by the respective diphoton branching fraction B[7]q 77]. For rjc, we may use the 
measured result BexpiVc — ^ 77] = 2.4+q;9 x 10""^ [ll]. For ric{2S), only an upper bound 
has been set experimentally, BcxpWc ~^ 77] < 5 x 10"^. However, there is good reason 
to believe B[ric{2S) 77] ~ B[ric 77]. To date the decay rjb — > 77 has not yet been 
observed experimentally. Theoretically, the branching fraction of rjQ — » 77 can be estimated 
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in NRQCD factorization approach 114"] 

r[^Q 77] 



2a2(2mQ) 1 + [(f - 5)C^ + (if - i|f^)C^ - |n,] : 

where we have approximated the total hadronic width of rjq by its gluonic width, and 
have included the NLO QCD corrections for both rig — > 77 and riq gg. Cp = 2N ' 
Ca = Nc, are the Casmirs for the fundamental and adjoint representations of SU{3)c group, 
respectively. Taking the number of active light flavors Uf = 4 for 77^, and as{2mb) = 0.18, 
we then get B['r]h 77] ~ 4.7 x 10"^. 

Owing to the energy conservation, rjc can only decay to e+e^ and fi^fJ^^, while rjd'^S) 
and rjb can access all three generations of leptons. From Table [H one can see that the values 
of r span a rather wide range, from the smallest 10~^ in r/f, —>■ e^e~ to the largest 0.96 
in 1]'^ T^T^ . Therefore, quarkonium decays to a lepton pair seem to provide a richer 
theoretical playground than the analogous decays of light pseudoscalar mesons. 

In the following we summarize the main lessons we have learnt from Table [H 

1. Among all the studied pseudoscalar quarkonium decays to lepton pair, r/f, t^t~ 
seems to have the largest branching ratio, ^ 5 x 10^^. However, the number of 
produced rjh at Super B experiment may not be copious enough for observing this 
decay mode. On the other hand, it also looks rather challenging to tag this decay 
mode in the high-energy hadron collider experiment such as LHC, since the r events 
are difficult to reconstruct in hadronic collision environment. It is also interesting to 
note that, the decay rjf, fJ'^f^~, with a branching ratio as small as 10~^°, seems 
comparable in cleanness with the decay chain rji, — > J/tpJ/ip — > 4/i, the "golden mode" 



for hunting 77^ at LHC, which has an estimated branching ratio of (0.7—6.7) x 10^^'^ |31 |. 
It may be worthwhile to look for this dimuon decay mode at LHC, but likely the signal 
events would be completely swallowed by the copious backgrounds. The rarest decay 
channels, rjb — > e+e~, 77c(?7c) — ^ e+e", with branching ratios about 10~^^ — 10"^^, seems 
completely out of the reach of any foreseeable experiments. 



Note that this formula works poorly for rjc decay. For a more satisfactory estimate of B[ric 77 ], it is 
important to resum a class of contributions associated with the running of the strong coupling as [30| . 
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2. For light pseudoscalar mesons decays to a lepton pair, one often resorts to the unitarity 
bound, which is obtained from cutting the intermediate photon hnes in the amphtude, 
to estimate the decay rate of P — /+/^. In some cases this simphfied but model- 
independent predictions seems not far below the exact results. Inspecting the phase 
pattern of /(r) in Table [H it is clear that this approximation can not make an accurate 
account for the majority of pseudoscalar quarkonia leptonic decay processes. 

3. The helicity suppression mechanism, which is manifested in the prefactor r in ([7j), 
plays a prominent role in dictating the size of each leptonic decay rate, and is much 
more important than the logarithmically running / function. This is clearly seen in 
the smallness of the R ratio for rjb — > e+e, even though the respective |/| is the largest 
among all the decay channels. In an alternative case, 7]c{2S) t^t~ , which hardly 
suffers from helicity suppression because of the rather large r, the respective branching 
ratio nevertheless remains small. This may be partly ascribable to the small |/|, and 
partly to the rather limited phase space available for this process, recalling the factor 
P contained in ([7j). 

4. The contribution of weak interaction is insignificant in rjclrj'^) decay, but can become 
important in rjiy leptonic decay. For example, including the Z° exchange effect will 
significantly enhance the branching fraction of 775 t^t~ predicted by QED alone, 
almost by one order of magnitude! This can be clearly understood from the fact that, 
the relative importance of the weak interaction effect grows with ttlq (see ([7])). The 
interference between QED and weak interaction can be either destructive or construc- 
tive, depending on the sign of the axial charges of heavy quarks, and also on the sign 
of the real part of the / function. 

Next we explore the phenomenological consequences of the Dalitz decay rjq —>■ l^l^'-f. In 
Table Hi] we tabulate the R ratios integrated over the full three-body phase space, as well 
as the corresponding R ratios by integrating over the photon energy up to 20 MeV. The 
purpose of including the latter is to assess the likelihood for these Dalitz events to fake the 
exclusive rjQ event. The main understanding we gained from Table HTl are: 

1. These three-body Dalitz decay processes in general have branching ratios several orders 
of magnitude larger than the corresponding rjq decay. This drastic disparity 



20 



TABLE II: The normalized R ratios for various pseudoscalar quarkonia Dalitz decays. R[E^ < 
20 MeV] represent the corresponding ratio with a 20 MeV cutoff imposed on the photon energy, 
so these Daltiz events may be experimentally indistinguishable from those r]Q events, 
i? [integrated] is the R value integrated over the whole three-body phase space. For comparison, we 
also juxtapose the SM predictions to the R values for pseudoscalar quarkonia decays to a leptonic 
pair, which are lifted from Table HI 



Decay modes 




r 








R[E^ < 20 MeV] 


ii [integrated] 


Tjc e' 


'e 7 




1.18 X 10" 


-7 


1.31 X 10" 


-9 


1.41 X lO"'^ 


0.0227 


Vc fJ'' 






5.03 X 10- 


-3 


1.82 X 10" 


-6 


1.41 X 10"'^ 


0.0062 


Vc{2S) - 




7 


7.90 X 10" 


-8 


9.63 X 10" 


10 


9.44 X 10"^ 


0.0233 


Vc{2S) - 




"7 


3.38 X 10" 


-3 


1.42 X 10" 


~6 


9.44 X 10-^ 


0.0068 


Vc{2S) - 


T^r" 


"7 


0.955 




6.35 X 10" 


-6 


2.73 X 10"^ 


2.8 X lO"'^ 


rjb e"^ 


'e~7 




1.19 X 10" 


-8 


3.07 X 10" 


10 


1.41 X 10"^ 


0.0263 


Vb 






5.07 X 10- 


-4 


1.29 X 10" 


-6 


1.41 X 10-^ 


0.0098 


% T^ 


^T~7 




0.143 




9.74 X 10" 


-5 


1.40 X 10-*^ 


0.0014 



should be attributed to several factors: less suppression by powers of a, the absence 
of helicity suppression, and the collinear enhancement. As a result, the branching 
ratios of f]c{Vc) ~^ e+e^7 may reach 4 x 10^^, and those of f]c{Vc) ~^ f^^f^~l 
reach 10~^. In the recently launched BESIIl experiment, roughly 10^'^ J/ip events are 
expected to be accumulated, rjc can be most copiously produced from J/ip through 
the magnetic dipole transition, with BcxplJ/'ip ~^ Vcl] = 1-3 ± 0.4% [11]. Therefore 
about 10^ rjc events are expected to be collected, and about (9(10^) Dalitz decay 
events should be produced. Even taking the detection acceptance and efficiency into 
account, the observation prospect for the aforementioned Dalitz decays at BES-111 still 
looks optimistic Similarly, if the future Super B factory can accumulate a huge 
rib samples, it may again be feasible to look for the rjb e"'"e^7, events. It is 

obvious that the majority of Dalitz decay events will obey a fragmentation pattern. 

However, the decay rj^ e+e^7 at BES-III experiment may be subject to substantial contamination from 
the background Bhabha events. In comparison, the decay rjc /Lt"'"/i^7 might be easier to establish. 
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i.e., a hard photon of ^ — ^ recoihng against a pair of nearly coUinear leptons. 

2. ?7c(25') T^T~'j is the only exceptional Dalitz decay process that has a even smaller 
decay rate than its respective leptonic decay rjci^S) t^t^. It is interesting to 
trace the reason. First by recalling f|T8|) . we note the rather strong suppression of the 
integrated decay ratio of the Dalitz decay near the mass threshold, oc (3^. By contrast, 
the suppression of the decay ric{2S) t^t^ is only linear in (3, stemming from the 
two-body phase space. It is these two very different threshold scaling behaviors that 
result in this anomalous pattern. 

3. Suppose a realistic electromagnetic calorimeter, say, the one installed in the BES-III 
detector, can detect those photons only with energy greater than 20 MeV |13|. When 
such a cutoff is imposed, according to Table [Tll the number of Dalitz decay events 
from ric{ri'J — > e"''e~7 with photon energy less than this cutoff, turns out to be about 
100 times greater than that from —>■ e^e~ . Therefore, what are experimentally 
recorded as the ric{i]'^) — > e~^e~ events, in fact receive the bulk of contribution from 
the three-body Dalitz decay events with unregistered soft photons. However, even 
if incorporating this two orders-of- magnitude enhancement, the corresponding decay 
probabilities are still too low for such channels to be established at BES-III experiment. 
A valuable lesson learned from this example is that, in any attempt to interpret the 
possible r]Q event as the signature for new interactions beyond SM, one must 
ensure that the contamination from the corresponding Dalitz decays has already been 
thoroughly understood and carefully incorporated in the analysis. 

V. SUMMARY AND OUTLOOK 

In this work, we have performed a comprehensive analysis of pseudoscalar quarkonium 
decays to a lepton pair without and with bremsstrahlung. These rare electromagnetic decay 
processes offer a clean platform to test our understanding of quarkonium dynamics. We 
corroborate the previous conclusion that the exclusive decays to lepton pair are extremely 
suppressed in Standard Model. By contrast, the pseudoscalar quarkonium Dalitz decays in 
general have a much larger decay rate, because of several joint factors: less suppression by 
power of a, absence of helicity suppression, and collinear enhancement. It is found that the 
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Dalitz decays r]c — ^ e^e~7 and rjc /i^/i~7, with branching fractions of order 10~^, may 
have the bright prospect to be estabhshed in the BES-III experiment. 

It is stressed that a reahstic electromagnetic calorimeter can detect photons only down 
to a minimum energy. Thus from the experimental perspective, those Dalitz decay events 
with photon energy less than this minimum energy, will be tagged as the exclusive lepton 
pair events. Taking this fact into consideration, it is found that the measured decay rate of 
Vciv'c) e"^e~ would be about two orders of magnitude greater than that literally predicted 
from ([7]). Nevertheless, such amplification is still not dramatic enough to warrant their 
observation at BES-III experiment. In general, the observation prospect for the pseudoscalar 
quarkonium decays to a lepton pair seems extremely pessimistic within the frame of the 
Standard Model. In this respect, future unambiguous sighting of any of this type of decays 
may be viewed as the strong evidence for the existence of new physics. 

Our analysis is based on a leading order calculation both in strong coupling and in 
quark velocity. To improve the reliability of our predictions, it is worthwhile to implement 
the QCD perturbative and relativistic corrections to the pseudoscalar quarkonium decay 
processes considered in this work, in particular to the quarkonium Dalitz decay processes. 

An interesting extension of this work is to study ^Pj {J = 0,1,2) quarkonium states 
decays to a lepton pair with and without bremsstrahlung. The processes Yca 2 — >■ e~^e~ 



and e~^e~ — > Xci,2 have been studied long ago within the color-singlet model 32j where 
infrared divergences are reported and subsequently cured by imposing a phenomenological 
cutoff of "binding energy" . To our knowledge, the process Xcj l^l^'y has not been con- 
sidered before, which is also plagued with infrared divergences. It is theoretically interesting 
to investigate these two processes in the context of NRQCD (NRQED), which provide a 
systematic way to tame these infrared divergences by incorporating the effect of higher Fock 
state (|cc(3^i)7) in our case, with the dynamical photon understood to be ultrasoft). We note 
that the latter process is quite analogous to the Xbj inclusive decays to charmed hadrons, 
which have recently been analyzed in the NRQCD factorization approach [33i]. Thus one 
may simply use their results with some slight modification. 

Eliminating the infrared divergences associated with the former process turns out to be 



For the calculations of Xci,2 ~^ e+e , one is allowed to put to since these processes are not subject 
to the helicity suppression, and retaining a nonzero mg only yields a small correction. 
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more subtle, and, more interesting, but likely to be achievable provided that one appeals 
to the even lower-energy effective theory of NRQCD, the so-called potential NRQCD (pN- 
RQCD) [34], by retaining the ultrasoft gluons (photons) as the manifest degree of freedom 
(for a pNRQCD-based study of the color-octet effect in exclusive reaction involving quarko- 
nium, such as J/tp ^ rj^j, see |37|). We hope future studies of these two types of processes 
at BES-III experiment, especially the resonant production process e^e^ Xci,2, may lend 
some important guidance. 
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